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1 Introduction 



There are not many exactly solvable potentials in quantum mechanics. As a result there have always been 
efforts to enlarge the class of exactly solvable potentials. Some of the different methods which have been 
used time and again to generate a hierarchy of isospectral potentials are the factorization method of Infcld 
and Hull , the Darboux algorithm , the method of supersymmetric quantum mechanics (SUSY QM) 
0, or the integral transformations of Abraham-Moses-Pursey etc. Among these methods the Darboux 
algorithm and the SUSYQM are closely related and these methods have found numerous applications in 
different areas of theoretical and mathematical physics 3 . 

At the same time, the scheme is still narrow as conventional SUSY fails to explain certain phenomena, 
e.g., the disappearance of the leading Borel singularity of the perturbation correction for the ground state 
energy of a SUSY theory 0. In order to explain such behaviour and also to widen the scope of SUSY 
QM, an idea was put forward to extend SUSY to higher orders 6 . We recall that in the conventional 
intertwining technique, two one-dimensional Schrodinger Hamiltonians H and H are intertwined by means 
of differential operators L as 

HL = LH HL^ = L^H (1) 

If L is of the first order in derivatives, the standard SUSY QM, with supercharges built of first order 
Darboux transformation operators, and the factorization method are recovered. On the other hand, if 
higher order differential operators are involved in the construction of L, it is variously referred to as 
polynomial SUSY 0, or nonlinear SUSY 0, or higher order SUSY (n-SUSY) 00, or Nfold SUSY 
0E3> the study of which has attracted the attention of a lot of researchers in recent times 00 El El 
EH . Contrary to standard SUSY, the anticommutator of the supercharges no longer coincides with 
the Hamiltonian in general. Instead, it becomes a polynomial of the Hamiltonian in degree N, and is 
sometimes referred to as the Mother Hamiltonian ll()j . 

Furthermore, the equivalence between an iV-th order Darboux transformation and a chain of N first 
order Darboux transformation is well established . Every chain of N first order Darboux transformation 
creates a chain of exactly solvable Hamiltonians ho — ► h\ /ijy. Hence the intertwining operator 

LS N ' between the initial Hamiltonian h$ and the final Hamiltonian can always be presented as a 
product of N first order Darboux transformation operators between every two juxtaposed Hamiltonians 
ho, hx,...,h,N : 

I w = L N L N - 1 ---L 2 L 1 h p L p = L p h p -i p = l,2,---N (2) 

In conventional higher order SUSY, ho and hjy are essentially self-adjoint Hermitian operators in a Hilbert 
space, with square integrable eigenfunctions. If all the intermediate potentials V\{x), ^2(2), . . . , Vn—i(x) 
are real valued functions in their common domain of definition (a, 6), the chain is called reducible, and 
the JV-th order Darboux transformation is called reducible as well. Additionally, if all the intermediate 
potentials are free of singularities in (a, 6), the chain and the corresponding transformation are called 
completely reducible. When at least one intermediate potential is a complex valued function, the chain 
and the corresponding transformation are called irreducible. 

At the same time, non Hermitian Hamiltonians have made an important place for themselves in the 
recent development of quantum mechanics, because of their intrinsic interest 1 1 1 1 and possible applications 
[T2*] . It is well known by now that a non Hermitian VT symmetric Hamiltonian admits real eigenvalues if 
the eigenfunctions, too, respect the VT invariance (the so-called unbroken VT symmetry), whereas the 
eigenvalues occur as complex conjugate pairs if VT symmetry is spontaneously broken (in this case the 
eigenfunctions are no longer VT invariant). For such non Hermitian VT symmetric Hamiltonians, 

VTH = HVT (3) 

where V stands for the space inversion operator and T denotes time reversal : 

V : x — > —x, p — > — p 
T : x — ► x, p — > — p, i — » —i 
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The reality of the spectrum may be attributed to the so-called ^-pseudo Hermiticity of the non Hermitian 
Hamiltonian ^3] 

fft = rjHr]- 1 (5) 

where n is a linear, invertible, Hermitian operator. Several non-Hcrmitian Hamiltonians, whether pos- 
sessing VT invariance or not, have been identified as r\ pseudo-Hermitian under rj = e~ 6p , where 9 is 
real, and p = — i-4-, or r) = e~^ x ', where 4>{x) is some gauge-like transformation. We note that for 
VT symmetric Hamiltonians, rj may simply taken as the parity operator V, whereas for conventional 
Hermitian Hamiltonians, rj = 1. 

Moreover, the square integrability of the wave functions is no longer a pre-requisite for non Hermitian 
Hamiltonians. Instead, the ortho-normalization of the wave function for Hermitian quantum mechanics 



^ V n dx = 5 m>n (6) 
is replaced by Jl] 

/ [CVT^ m ]^ n dx = 5 (7) 
where C plays the role of a linear charge operator, obeying the relationship 

[C,H] = [C,VT]=0 (8) 
and has the property C 2 = 1. In the position representation C is given as 

C(x,y) =22ip n (x)ij} n (y) (9) 

n 

and the completeness relation gets modified to 

J2 [CVT^ n {x)\ Mv) - S{x - y) (10) 

n 

While nonlinear SUSY for N = 2, has been investigated widely for Hermitian Hamiltonians [5-10], such 
studies have not been carried out as yet for non Hermitian Hamiltonians. Motivated by the importance 
of such systems in the recent development of quantum mechanics, our aim in the present work is to 
generalise the concept of nonlinear SUSY to include non Hermitian quantum systems. In analogy with 
the first order systems, where the partner Hamiltonians H± of non Hermitian systems were found to be 
related through pseudo super symmetry [181 115) , it will be shown that the underlying symmetry between 
the isospectral partners ho and hjsr is a generalisation of TV SUSY and may be called nonlinear pseudo 
super symmetry. The nature of the intermediate Hamiltonians as well as the corresponding wave functions 
will also be investigated. 

The organization of the paper is as follows. For the sake of completeness, in section 2 we briefly 
outline conventional nonlinear SUSY for Hermitian quantum mechanics. In section 3 we describe a 
similar framework for non-Hermitian Hamiltonians and show that, the underlying symmetry for the 
potentials produced by higher order Darboux algorithm, is nonlinear pseudo supersymmetry. Some 
explicit examples are given in sections 4 and 5 , while section 6 is devoted to a conclusion. 

2 Non linear SUSY for Hermitian Hamiltonians 

In the conventional first order supersymmetric quantum mechanics, if a given solvable Hamiltonian 

H—^ + VM (11) 
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possesses a discrete spectrum of bound states E n , n = 0, 1, 2, • ■ •, together with the square-integrable 
eigenf unctions ip n (x) , then a pair of first-order operators Lq and Lq can be constructed from the ground 
state tpo, given by 



Lo = Tx+ W ( X ) 



JQ ~-dx~ + W ° ix) (12) 
where 

W (x) = -[ln^ (x)]' (13) 

such that Lq and Lq play the role of intertwining operators for the initial and final Hamiltonians H and 
H, respectively : 

HL = L H HLl = l\h (14) 

with 

H = l\l , H = L Ll (15) 
Simple straightforward algebra shows that the partner potentials V(x) and V{x) can be expressed as 

V(x) = W*(x)-WZ(x) 



V(x) = W${x) + W^x) = V(x) + 2W^(x) 
The eigenfunctions ^{x) and ip(x) of H and H are interrelated through Lq and ij : 



L ^ Q {x) = 0, L (x)ipi(x) = W "^ x ^ a ipi^x), 4i ) i( x ) a ipi(x), 

IpoyX) 



i = l,2,- 



(16) 
(17) 

(18) 



where Wo,i(x) — {ipo(x)ip'i(x) — tp' {x)il)i{x)} is the Wronskian of ipo(x) and ipi(x). The concise algebraic 
form of spectral equivalence is given by the superalgebra for the partners H and H , and the supercharges 
Q and : 

' L 





H = {Q,Q^} = 



f = 





4 o 



( H 


> 




\ o 


H j 


■( 



Lq4 



satisfying the relations 



{Q,Q} = {Q t ,Q t } = 0, [Q,H]=[Q\H]=Q 



(19) 
(20) 

(21) 



Thus H and H are isospectral except for the lowest eigenvalue Eo which is missing in H, as tpo is not 
normalizable. 

To generalize standard SUSY to higher order, the supercharges are built of higher order intertwining 
operators [S]. The two Hamiltonians ho and ft. at are intertwined through an N-th order differential 
operator as 

LWfto = %LW, ft i (W)t =i (W)t ^ (22) 

where /io and hjq are self adjoint operators. The proper eigenfunctions tpi of the original Hamiltonian 
are known exactly : fiQipi = Eiipi. Any such operator can always be presented in the form known 
as Crum-Krein formula |16| 



(23) 
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=W- 1 (u 1 ,U2,■ 
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where W (u±, 112, ■ • ■ , ujv) stands for the usual symbol for the Wronskian of the functions u\, 112, ■ ■ ■ , un- 
The functions Uj, (z = 1,2,..., iV) called the transformation functions are eigenfunctions of /iq, /low, = 
OLiUi, and they need not necessarily satisfy any physical boundary condition. The final potential has the 
form 

V N (x) =V(x)-2-^]nW( Ul ,u 2 ,...,u N ) (24) 

and will be free of singularities whenever the Wronskian is nodeless, which in turn, requires that only 
consecutive eigenfunctions of ho have to be considered The eigenfunctions ip{x) and ipi x ) °f ^0 and 
/ijv are connected by the intertwiners L^ N > and L'"'^ as : 



W "j,j+l,-,j+iV,i( a: ) 

^j,j+i,-,i+Jv(a;) 



(25) 



where VV}j+i.....j+jv,t(aO and VFjj+i.....j+Ar(a;) are the Wronskians of the eigenfunctions of ho associated 
with the corresponding subindices. Thus if ipi(x) is an eigenfunction of ho with energy E i: then ipi(x) is 
an eigenfunction of h^ with the same energy E^. Evidently 



I ( % = 0, i=l,2„---,iV 
However, for energies = 1, 2, • • • N), the corresponding eigenfunctions of hjq 

ip(x) 



(26) 



ip(x) a 



Wj,j+i,...,j +N (x) 



have growing asymptotics at both infinities. Consequently these are not physically acceptable solutions 
of ft.jv, and the corresponding eigenvalues Ei(i = 1, 2, ■ ■ • N) are excluded from the spectrum of hN- Thus 



(27) 



with the exception of the levels E = Ei, i = 1, 2, • • • , N, which will be absent in the spectrum of the new 
Hamiltonian hx, as the corresponding eigenfunctions are not square integrable. 

It has already been shown J7| that the operator l'"' can always be presented as a product of N first 
order Darboux transformation operators between every two Hamiltonians ho, h\, . . . , h^ : 



= LnLn—i ■ ■ ■ L\ hpLp = Lph 



pilp—X 



p=l,2,..-,JV 



(28) 



We note that the final Hamiltonian hjsr is Hermitian, although some of the intermediate Hamiltonians 
hi could be unphysical, e.g., their associated potentials might contain extra singularities that were not 
present in the initial one. The supercharges Qn and Q N are constructed as 



Q 



N 



LW 




Ql 



N 





L Wt 



Evidently, Qjv and Q N are nilpotent 



{Qn,Qn} = {q n ,Q n } = 



The super Hamiltonian 



satisfies the relations 



H 



N 



ho 






h N 



[Qn, Hn} — Qjy, H 



IN 



= 



(29) 

(30) 

(31) 
(32) 
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The anticommutator can be generally expressed by a iV-th order polynomial Vtf of the Hamiltonian Hn 

Km = jQ^Qjv} = ( L q l {n) l (n)\ ) = n ( Hn ~ akI ^ ( 33 ) 
^ ' fe=i 

where X is the 2x2 unit matrix, and 

N 

LWt£W = J] (/l0 -a fe ) (34) 

k=l 
N 

LWL^=l[(h N -a k ) (35) 
fe=i 

Since the right hand side of l|33|l is a polynomial in H^, it is called nonlinear SUSY or iV-fold SUSY. 
The operator 7ij\f is termed as the Mother Hamiltonian and satisfies the commutation relations |U] 



[Qn, 7~i 



Q^ N ,H N ]=0 (36) 



For N = 1, iV-fold SUSY reduces to standard SUSY. 

The most widely studied higher order SUSY is for N = 2 [HUH!) where the formalism reduces to : 

i (2) = L 2 L X (37) 

where 

L x = -d x + Qjiu 1 )', L 2 = -d x + Q.TL V y, v = L lU2 (38) 
and the isospectral potential turns out to be 

V 2 (x) = V(x) - 2^ In W jtj+1 (x) (39) 



3 Non linear pseudo-SUSY for non Hermitian Hamiltonians 

In this section we extend the concept of nonlinear or iV-fold supcrsymmctry to non Hermitian quan- 
tum mechanics. Though the Darboux algorithm and (nonlinear) supersymmetric quantum mechanics 
are equivalent for Hermitian Hamiltonians, the situation is different for non Hermitian Hamiltonians. 
However, intertwining operators and B^ N ' can still be constructed with the help of Darboux trans- 
formation. Analogous to the case of Hermitian quantum mechanics, it will be shown that once a non 
Hermitian Schrodinger potential V(x) is exactly solvable, one can construct an isospectral partner Vn(x) 
from (1231) 

V N (x) =V{x)-2-^hxW( Ul ,u 2 ,...,u N ) (40) 

where W stands for the usual symbol for the Wronskian of the functions u±, u 2 , . . . , itjv, which are 
eigenfunctions of ho, hom = a%Ui. As before the functions Ui(x) may be just formal eigenfunctions. 
Our aim will be to study the spectrum of the new Hamiltonian in detail, to investigate the nature of the 
potential and the eigenfunctions, and to determine the symmetry which connects the original Hamiltonian 
ho and the transformed one hjy. For this purpose, we look for two intertwining operators A^ N ' and B^ N ' 
such that 

AWfco = M W , h BW=B^h N (41) 
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where ho and hw are no longer self-adjoint operators (/io,(jv) rm) 5 on the contrary, to ensure the 
reality of the spectrum, they are rj pseudo Hermitian 

vho.(N)^ 1 = h f 0{N) (42) 

where r) is a linear, invertible, Hermitian operator. However, the choice of r\ is not unique. For VT 
invariant potentials, a simple representation of rj may be given by the parity operator : 

r, = V, Vf(x) = f(-x) (43) 

It follows that for real potentials, Ij43() leads to rj = 1 so that = A^ N '^, thus reproducing the standard 
result of supersymmetry. 

It follows from equations (|41|l and l|42|) that the operators and B( N > are pseudo-adjoint : 

B (N) = A (N)# = V -1 A (N» V (44) 

Considering first order Darboux transformation between every two juxtaposed Hamiltonians ho, hi, ... , hjy, 
each pair intertwined by first order operators L k (k = 1, 2, • • • , N) 

h k L k = L k h k _ x fe = l,2,3,---,JV (45) 

L*h k = h k ^L* k = 1,2,3,- ~,N (46) 



if = 77-^^ (47) 



where 

\ 

then, analogous to the Hermitian case, the final Hamiltonian hjsr is found to be related to the initial (or 
starting) Hamiltonian ho through 

h N = L N L N -i---L 2 L 1 h a L*L*---L* (48) 

so that the operator A^ N ' can be represented as a product of the N first order Darboux transformations 

AW =L n L N -i...L 2 Li (49) 

with its pseudo-adjoint 

gW = A{ N)# = Ln7] = L # L #_ iL # (5Q) 

It is worth mentioning here that in contrast to Hermitian quantum mechanics, all the intermediate Hamil- 
tonians h k are physically acceptable as their associated potentials contain no extra singularities which 
are not present in the initial potential V(x). This is essentially because the associated eigenfunctions do 
not have nodes on the real line, and they are normalizable in the sense of equation l|7)l. 

Thus the initial and the transformed Hamiltonians ho and hx are related by non linear pseudo su- 
persymmetry. The super Hamiltonian of this system consists of the pseudo supersymmetric pair of 
Hamiltonians ho and hN as 

Hn = ( H ° °) (51) 



h 7v 

The supercharges generating this form of pseudo-supersymmetry are constructed in the following way : 
AW A ^ _i_ t { W 



n = I J , Qjj = V Q J N V = I BiN ) o ) = I A (N)# o J (52) 
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so that the supercharge Qn and its adjoint Q N of standard Hermitian quantum mechanics are replaced 
by Qn and its pseudo adjoint Q N f° r non Hermitian Hamiltonians. Obviously, Qn and are nilpotent 

{Qn,Qn} = {q%,Q%}=0 (53) 

and satisfy the following closed algebra : 

[Qn, Hn] = [Q%,H N ] =0 (54) 

r 4, ~\ ( \ A 

W;v = |Q#, Qat} = = II - (55) 

' fe=i 



i.e., 

AT 

AW#AW = n(fto-a fc ) (56) 
fc=i 

N 

AWAW* = f[(h N -a k ) (57) 
fe=i 

and T is the 2x2 unit matrix. Evidently, if i\)%{x) is an eigenfunction of fto with energy eigenvalue Ei, 
then ^j(a;) = A^ip^x) is an eigenfunction of Hn with the same energy Ei. However, for i = 1, 2, • • • , N, 

Hx) a (58) 

w{i} 1 ,i} 2 ,---,^N) 

Clearly, the eigenfunctions tpi(x) (i = 1,2,- ■ ■ N) of ft, at corresponding to the eigenvalues Ei (i — 
1, 2, • • • , TV) grow asymptotically, and so cannot be included in the set of solutions of ft. at. Consequently, 
Ei (i = 1, 2, • • • , N) are excluded from the spectrum oi h^- 

Next we note two interesting results which are in contrast to the Hermitian case: 

1. For Vn (x) to be free of singularities, the Wronskian W (^1,^2, V'jv) = W^ 1 ^ 2v ...^ N (x) must 
be nodeless. In case of Hermitian potentials, this is guaranteed only when tpi, i — 1, 2, • • • N represent 
N consecutive eigenfunctions. However, in case of generic non Hermitian potentials, the eigenfunctions 
i> n (x), (n = 0, 1,2, •••) have no nodes on the real line. Consequently, the Wronskian is free of real 
singularities for any value of i,j,k,---, and thus can be used to generate a wider class of isospectral 
Hamiltonians. 

2. The intermediate Hamiltonians are also physically acceptable, as the corresponding potentials are 
free of singularities, for the same reason as given above. For example the first intertwining gives 

V 1 {x) = V{x)-2^\ni, i {x) (59) 

which is well defined. However, this may not always be true for Hermitian potentials due to the presence 
of additional singularities in V\{x), which are not present in V(x). 

For the sake of simplicity, in the present work we shall restrict ourselves to second order nonlinear 
pseudo-supersymmetry. Thus if an intertwining operator A = L^Li is constructed from the two first 
order Darboux transformation operators L\ and L2, given by 

Li = -d x + (lnitj)', L 2 = —d x + (lnu)', v = L^Uj (60) 
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where and Uj are any two eigenfunctions of the non Hermitian Hamiltonian ho, then the transformed 
isospectral Hamiltonian 

h2 = -^ + Vij(x) (61) 

has eigenfunctions 



= -E n ip n + E l i) l — — — f- + Ejipj — — — — 

where 



d? 
dx 

The mother Hamiltonian Hi is constructed from the anticommutator by 



(62) 



Vij(x) =V(x)-2-^lnW(ui,Uj) (63) 



*. - {Qt.Q,} - ( a : a & ^ - - a - i] , h - *,x>V - 02 x, ) <«> 

where X is 2 x 2 unit matrix and ffa is given by l|51|) . 

In the following sections we shall investigate this formalism further with the help of explicit examples. 



4 VT Symmetric Oscillator 

In this section we shall apply our formalism to the well known example of the VT symmetric oscillator 

M 

2 1 

V(x) = (x- it) 2 + a ~* (65) 
(x — te) 

with eigenfunctions 

i> n (x) = e -*<'- fc > a (x - ie)- qa+ ^ L- qa ((x - ie) 2 ) (66) 

and eigenvalues 

E n = An - 2qa + 2, n = 0,1,2,- •• (67) 

where q = ±1 is called the quasi-parity. 

In this study we shall restrict ourselves to N = 2 only. If one performs Darboux transformations with 
two eigenfunctions ipi(x) and ipj{x) of the potential V(x), corresponding to energies Ei and Ej (i and j 
need not be consecutive), then the intertwining operators take the form 

1 dx tpi ' 1 dx ipi ^ ^ 

l d M t*- d + w * ^ (m) 

where Wij is the usual Wronskian given by 

Wij = Wtyi,^) = ^i{xWAx) - ^(x)^(x) (70) 
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and r\ has been taken as in (|43|l . Replacing the intertwining operators A^> (and B^) by A (and B) for 
simplicity, we obtain 

A-LM- ^±-^&-& (71) 

'-^ - ,72 > 

where 

W' ■ 

and 77 is simply the parity operator V for VT symmetric potentials. The new exactly solvable non 
Hcrmitian potential, which is isospectral to the VT symmetric oscillator in Ij65(l . is obtained from 

V id ( x ) = V{x)-2^\nW id (74) 

with solutions 

$ k {x) = A?p k (x) (75) 

Thus for each set (i,j), one obtains two sets of Vij (x) because of the presence of quasi parity q. Obviously, 
ipk(x) = for k — Thus the new potential so constructed, in l|74|) above, has all the eigenenergies 
of the original VT symmetric oscillator except for the levels i, j, which are missing from the spectrum of 

m- 

For the simplicity of calculations we shall now construct and examine some potentials using low values 
of i and j in further detail. 



4.1 New potential for i — 1, j — 2 

Applying the above formalism with the two eigenstates ipi(x) and ^ 2 (x), of the potential in (|65|) . the 
Wronskian is found to be 

W (Vi, V2) = c 12 e-(— ' 2 (x - ief- 2qa g (76) 
where C12 is some real constant and 

g = (1 - qa) (2 -qa) -2(1- qa) [x - it) 2 + (x - ie) 4 (77) 

The intertwining operators A and are obtained from A = L2L1, — LfLf, where 



Li 



dx Vl f7ftl 
d , , -qa+ 1 2(a;-ie) ^ 

— ; (a; - ie) + r-f 9 

ax (x -ie) i-q a -(x- i e ) 

d d . , s 

d , , -ga+| 2(x-ie) 0' ( 79 ) 

•— - (x - ie) + -7 r~r" + 5 + " 

(a; - ie) 1 _ qa - (x - ie) 9 
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so that 



.4 . ^ H -,(,- fc) + ^ + gu + ( ,- fc) v ^t^T + >) 



da; 2 [ (a; - ie) g J da; ( x _ i e ) 2 

-get +5 2 (a; — ie) 1 



(80) 



+ 2qa -1+ i-(x-ie) 



(x ~ie) 1 - qa - (x - ie) \ 9 



A* = jl + {-2 (x - ie) + 2 g -^ + *}l + ( X - ief + ^ + t ® ;T " *> 



cb 2 [ (x - ie) g j dx ( x - ie) 

o f / x -9" + ^ 2(x-ie) 1 g' g" /#' x 2 
+ 2ga-3+<J - (x - ie) + — j-f - : i — . 2 ' 



(81) 



(2 - «0 1 - get - (x - ie) J 5 5 
Applying equation l|74|l. the new potential isospectral to the one in Ij65|) except for the states corre- 



sponding to ip\{x) and ip 2 {x), comes out as 



Vi, a (a) = (x - ief + K - 1 \ 2JK \~ ' 2J - 2— + 2 I — ) +4 (82) 

(x - ie) 5 V 5 / 

which has solutions 

T / s P , , p , W 7 ('0„+2,V'2) , W(^l,V'n+2) / S oN 

with energy eigenvalues 

E n = E n+2 = 4n + 10-2ga , n = 1,2,3, ■•■ (84) 

The ground state is given by 

^o(x) = e-H*-*>\x ~ ie)-^+* |s x + ^^"^ j (85) 

with eigenvalue 

E = E Q = 2- 2qa (86) 

where -Bi and -B2 are some x, e independent constants. Thus the energies E\ and -E2 of F(x) are absent 
in the spectrum of V\ l2 (x). It can be verified that the eigenfunctions tp are also VT invariant, and can be 
normalized using Q. Furthermore, the supercharges Q 2 and Q 2 , generated from the operators A and 
A#, satisfy the following algebra : 

H 2 = {Q2, Qf } = Hi - 4(4 - qa)H 2 + (6 - 2?a)(10 - 2ga) (87) 

where H 2 is given by (|51|l . The intermediate potential given by 

Vi{x) = V{x)-2^hx^{x) 

= (x-ie) 2 + t^±iHz^±l) + ^ ^ 

(x-ie) 1 - $a - (a; - te) |l - ?a - (. T - ze ) 2 | 

(88) 



(-ga+|) (-ga + |) o" , /g 



2 
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does not have any singularity on the real line, and hence is physically acceptable as well. By arguments 
similar to those given above, its ground state eigenfunction is given by : 

' e-^ x - te)2 (x-ie)- qa+ i (89) 



1 — qa — (x — ie) 2 
with energy 

e = E Q = 2 - 2qa (90) 

and the excited states 

<Pn = ; 91) 

01 

with corresponding energies 

e„ = E n+ i = An + 6 - 2qa, n= 1,2,3, ■■■ (92) 

It is easy to observe that applying (J3J, both the intermediate and the final potentials (as well as their 
eigenfunctions) satisfy ©, and hence are VT invariant, having real spectra. 



4.2 New Potentials for % = 0, j = 2 



In a similar manner, the expressions for the different quantities are obtained as follows : 

W (0 O , 2 ) = coae-<- fc > a (x - ie)~ w { %^ - 1 ) 

\ 2-qa 



with C02 some real constant 



Li = -— - (x - ze) + ±- rf- 

dx (x — ie) 



L 2 = -~r - (x - ie) + 



2(x-ie) 



A 



dx 2 



dx 



+ 2{(x-ie) + 



(x — ie) 
(qa - 1) 



(x — ie) — (2 — qa) 
(x — ie) I d 



+ 



(x - ie) (x - ie) 2 - (2 - qa) J dx 

(-«« + !) (-«« + §) 3 



+ (x - ie) z 



(x — ie) 2 



(x — ie) 2 — (2 — qa) 



2qa - 3 



(x — ie) 



A* = ^ 2 {(x -c) I > — 

dx 2 ) (x - ie) (x _ i e f - (2 - qa) J dx 



(x-ie) 2 {(x-ie) 2 - (2-qa)} 



4(-qa + 2) 



{(x-ze) 2 -(2- g a)}' 



2qa - 5 



(93) 

(94) 
(95) 

(96) 



(97) 
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The new potential 



V ,a(x) = (x-teY + T ^ — -f + , 



(a; — ie) (x — ie) — (2 — ga) 
8(2-,*) ^ 



{{x-ief -(Ji-qa)}' 

where 



5 

a = -qa+- (99) 

is totally different from the initial potential of the VT symmetric oscillator, yet shares the same spectrum 
except for the states n = 0, 2 of the original potential, which are missing in the partner. 
The ground state wave function of the Hamiltonian in (|98|l is given by 

Mx) = {A 1 (x- lef + A 2 + -A- 1 e -s(*-<« ) 2 (x - ie)- qa+ i (100) 

[ (x - ie) - (2 - ga) J 

with ground state energy 

E = Ei = 6 - 2qa (101) 
where A\, A 2 , A3 are ^-independent constants, while the excited states are obtained from 1)62(1 

rp n = A^p n+2 

^(-00,-02) ^(^0,02) 

with energies 

j?„ = £ n+2 = 4n + 10 - 2qa , n=l,2,--- (103) 

It can also be verified that eigenfunctions ijj n (x) have correct asymptotic behaviour and are also VT 
invariant. Consequently, they also satisfy equation J7J. The intermediate potential is given by 

V l{ x) = (x - *e) 2 + + ^ (T + l) + 2 (104) 

(x — ie) 

which is also physically acceptable. By arguments similar to those given above, its ground state eigen- 
function is given by 

O = e -*(*-«) a ( x - te )- qa+ i (105) 

with energy 

e = 6 - 2qa (106) 

and excited states 

W (ip n+u ip ) 

9n = ; (107) 

ipo 

with energies 

e„ = E n+1 = 4n + 6 - 2qa , n = 1,2,3, ■•• (108) 

Once again, both the intermediate and the final potentials (as well as their eigenfunctions) are VT 
invariant, having real spectra. 
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The supercharges Q2 and Q\ generated from the intertwining operators A and can be shown to satisfy 
the following algebra : 

n 2 = {Q 2 ,Q*}=Hl-4(3-qa)H 2 + (2-qa)(10-2qa) (109) 

where H2 is given by (|51|l . 

We note that the potentials obtained in this section are unique in the sense that they do not have any 
counterpart in standard quantum mechanics (i.e., in the Hermitian case). . 



5 VT symmetric Scarf II potential 

We note that the generalised oscillator problem considered in the last section was made non Hermitian by 
an imaginary displacement of the coordinate variable x. However, there are other methods of constructing 
non Hermitian models. To see how the formalism described in section 3 works with such models, in this 
section we shall study an example, viz., the VT symmetric non Hermitian Scarf II potential, which has 
been VT symmetrized in a different way. This exactly solvable potential, given by 

V(x) — —A sech 2 x — i/i sech x tanh x, A > 0, [i ^ (HO) 

has a discrete spectrum that admits both real as well as complex conjugate energies, depending on the 
relative strengths of its parameters A and fx. For < A + j, the system possesses a real and discrete 
bound state spectrum, whereas for > A+ |, the system exhibits spontaneous VT symmetry breaking, 
with complex conjugate pairs of energies. The normalized wave functions for this potential are well 
known, being given by |151 ITU] 

nW ( 5 - 2p) 

where are the Jacobi polynomials |2U| : 

P^(i sink x) = -E^±^±iL- F (-n, n + a + + 1; a + 1; z) (112) 
1 (n + 1)1 (a + 1) 



and 



1 — i sinh 



2 










h A - 


f /i = 








l\- 


hA 


-M = 



(113) 



Q =-i ± iV2 +A -^=-i ± l (115) 

However, for normalization of the wave functions, only the positive sign is allowed in p. The energy 
spectrum 



E n = -{n-p-qf , 7t = 0,l,2,... < ^ ' 2 J (116) 

is real and bound for < A + j , i.e., for real p and q, with two towers characterized by the two values 
of q. 



s + t-l 
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If the formalism developed above is applied to this example for N = 2, with states ?po(x) and ?p2(x), 
then the Wronskian is calculated to be 

W (ipQ, 1P2) = (1 — * sinhx) _2p (1 + i sinh x)~ 29 cosh 2; |— i(p — q) + ^p + q — sinhxj (H7) 

and the intertwining operators A and are given by 

A = L 2 L X A* = L*L* (118) 



where L\ and L2 take the form 

Li 



dx tpo 

— h i(p — q)sech x — (p + q)tanh x 

dx 



(119) 



L 2 = d 1 W ^ 2 < 
dx W 0t2 ipo 

— — — h i(p — q)sech x — (p + q) tanh x (120) 

(— p — q+ |) + i(p — q) sinh x + (—2p — 2q + 3) sinh 2 a; 
i(p ~ q) cosh £ + (— p — q + |) sinh a; cosh x 

Now using (|40(l the new potential is found to be 

~ ~ 2 „ / er 2 sech 2 x — i p a sech x tanh x\ . . 

m),2(^) = — A seen, x — 1 fi seen x tanh x — 2 = (121) 

Y (p seen, x — i a tanh x) J 

where 

A = A-4p-4g + 2 (122) 

/l = /i-4p + 4g (123) 

A = 2(p 2 + (? 2 ) + (p + (? ) (124) 

/! = 2 (p 2 - g 2 ) + (p - qr) (125) 

p = p-g (126) 

-P-9+I) (127) 

Once again, the final potential Vq$(x) is also VT invariant. The eigen functions are obtained from I|tj2|) . 
with the ground state as 

i> = {E Q -E 1 )^ 1 + (E 2 -E Q )^ 2 ^ (128) 

and excited states ^ 

i> n = (E - E n+2 ) ^ n +2 + (E 2 - E ) ^2%2 (129) 
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— 2p— — —2q— — 

where P n denotes the Jacobi polynomial P n 2 ' 2 (i sink x) and P' n denotes its derivative with 
respect to a;. It can be shown that for < A + j, the wavefunctions ip are also VT invariant, and can 
be normalized following (JJJ). The new potential Vo^(x) has real bound state spectrum given by 

E Q = -(l-p-q) 2 (130) 



E n = -(n + 2-p-q) 2 , n = l,2,---,< ( J ' "- - j (131) 

and the algebra satisfied by the supercharges turns out to be 

H 2 = {Q2, Qf} = H 2 + (2-2p- 2q) H 2 + (p + q){p + q-2) (132) 

where H 2 is given by 1)510. 

The intermediate potential takes the form 

V\{x) = —v\ sech 2 x — iv 2 sech x tanh x (133) 



s + t-5 



where 



with eigenfunctions 



v 1 =X-2(p + q) (134) 
v 2 = n - 2 (p - q) (135) 



(136) 



and the corresponding energies 

e„ =E n+1 = -{n + l-p-qf , n = 0, 1, ■■■,<( "' " - I (137) 



Thus Vi(x) and the corresponding wave functions l|136fl are also physically acceptable as well as VT 
invariant. 



6 Conclusions 

In this article we have suggested an application of higher order Darboux algorithm to non Hermitian 
VT symmetric potentials. For the sake of defmiteness the method has been applied to two specific 
potentials, namely, the generalised oscillator and the Scarf II potentials and a number of new potentials 
having nearly the same spectrum as the original ones have been obtained. It may be noted that in 
each of these cases, starting from a VT symmetric potential we have obtained new potentials which are 
again VT symmetric. In other words the higher order Darboux algorithm does not induce spontaneous 
VT symmetry breaking. Among the different cases considered here the one involving non consecutive 
levels deserves special mention. The potentials thus obtained have no Hermitian analogues. Also the 
intermediate potentials in all the cases are perfectly well behaved since the Darboux algorithm does not 
introduce any new singularity or break VT symmetry. Furthermore it has been shown that the symmetry 
underlying the original and the new potentials is a fusion of nonlinear SUSY and VT symmetry which 
we call nonlinear pseudo super symmetry. Finally we note that analogous to the study of breaking N fold 
supersymmetry |21) . it would be of interest to examine breaking of this new symmetry. 
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